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Abstract. We introduce a new method for computing massless Feynman in- 
tegrals analytically in parametric form. An analysis of the method yields a 
criterion for a primitive Feynman graph G to evaluate to multiple zeta val- 
ues. The criterion depends only on the topology of G, and can be checked 
algorithmically. As a corollary, we reprove the result, due to Bierenbaum and 
Weinzierl, that the massless 2-loop 2-point function is expressible in terms of 
multiple zeta values, and generalize this to the 3, 4, and 5-loop cases. We find 
that the coefficients in the Taylor expansion of planar graphs in this range eval- 
uate to multiple zeta values, but the non-planar graphs with crossing number 
1 may evaluate to multiple sums with 6 th roots of unity. Our method fails for 
the five loop graphs with crossing number 2 obtained by breaking open the 
bipartite graph Kz,A at one edge. 



1. Introduction 

Let n\, . . . ,n r 6 N and suppose that n r > 2. The multiple zeta value is the real 
number defined by the convergent nested sum: 

1 



C(ni,...,7i r ) = ^2 



0<k 1 <k 2 <...<k r 1 

An important question in perturbative quantum field theory is whether multiple 
zeta values, or some larger set of periods, suffice to evaluate a given class of Feyn- 
man integrals. Moreover, it is crucial for applications to find efficient methods for 
evaluating such Feynman integrals analytically. 

In this paper, we shall consider massless Feynman integrals in </>| theory with 
the propagators raised to arbitrary powers, and a single non-zero momentum. The 
simplest case is the two-loop example, pictured in figure 1 on the left. Let be 
a positive real number corresponding to each edge 1 < i < 5. The corresponding 
Feynman integral is: 

where D is the number of dimensions, and q is the momentum entering on the 
left. Now suppose that a t = 1 + n,- s, where rij are positive integers, for 1 < 
i < 5, and e is the parameter in dimensional regularization, i.e., D = 4 — 2e. 
The problem of calculating the coefficients in the Taylor expansion of |T]) with 
respect to e, is important in three and four loop calculations, and has a history 
spanning approximately twenty- five years (to which we refer to [3] for an account). 
In particular, it had been conjectured for a long time that every coefficient is a 
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rational linear combination of multiple zeta values. This question was finally settled 
in the affirmative in 3J, using Mellin-Barnes techniques. 

Until now, however, there seemed to be a lack of systematic methods for comput- 
ing a range of Feynman integrals analytically at higher loop orders. In this paper, 
we introduce a new method, using iterated integration with poly logarithms, which 
was initiated in [8] to compute the periods of moduli spaces of curves of genus 0. 
Using this, we reprove the fact that the Taylor expansion of the massless two-loop 
two-point integral (p} evaluates to multiple zeta values, and extend the result to 
higher loop orders. Our method also yields results for certain examples of massive 
Feynman diagrams, but in the present paper we only consider massless cases. 

1.1. Results. We consider three and higher-loop integrals with exactly one non- 
zero momentum and arbitrary powers of the propagators, generalizing the integral 
(JTJ) . Let G be a graph in <j>\ with two external legs, and a single momentum q. 
By Hopf algebra arguments, we can assume that G is primitive, in the following 
sense. Let G denote the graph formed by closing the two external legs of G, which 
now has no external edges, but gains an extra loop (see figure 1). We will suppose 
that G is primitive divergent in the sense of [4], i.e., G contains no strict divergent 
subgraphs and satisfies: 

#edges of G — 2 x #loops of G . 

In this case, we say that G is broken primitive divergent (bpd), and every bpd 
graph with a fixed number of loops can be obtained by taking the set of all primitive 
divergent graphs with one more loop, and breaking them open along every edge. 

Now suppose that G is bpd with h loops and L internal edges. It follows from 
the above that L = 2h + 1. For each internal edge i, let aj = 1 + n.i e, where is 
a positive integer, and D = 4 — 2 e, and consider the massless Feynman integral 

(2) /rt.../^ L n(i) Q \ 

2—1 * 

where ki is a momentum flowing through the i th edge, and is the propagator cor- 
responding to the i th edge. The domain of integration is given by the conservation 
of momentum at each vertex, and there is a single external momentum entering G, 
denoted q. It is easy to show, using the Schwinger trick, that (|2|) is proportional to 
a certain power of q 2 pj] . The constant of proportionality is given by a product of 
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explicit gamma factors and powers of n, with the parametric integral: 



(3) 



1(G) = f H 




doei . . . da\ . . . dotL+i , 



where Ug(at, . . . , Ol+i) is the graph polynomial (or Kirchoff polynomial) of G, 
at, . . . , Q-l+1 are Schwinger parameters for the L + 1 edges in G, A is any index 
between 1 and L + 1, and is a certain linear combination of at, ■ ■ ■ , cll- The 

problem of computing the Taylor expansion of ((2]) reduces to computing the expan- 
sion of ((3]). For example, if G is the wheel with 3 spokes, we have: 

— aia2a6 + aia4a6 + a2Ct5a6 + a4a5a6 + a i a 3 a 6 + a 2^3^6 + a 3 a 4 a 6 + Q; 3 Q; 5 Q; 6 
+ aia3a4 + Q!ia3Q!5 + a2a3a4 + a2a3a5 + a2Q;4a5 + aiQ;4a5 + aia2a5 + Q;ia2a4 , 

and the coefficients in the Taylor expansion of ^ are given by the period integrals: 



where mi, . . . , itil+i, n are arbitrary positive integers. 

We know state some results on the transcendental nature of the coefficients in 
the Taylor expansion of 1(G) with respect to e, for all bpd graphs G up to five 
loops. For every graph G for which a theorem is stated below, there is also a 
corresponding algorithm for computing the coefficients of the Taylor expansion of 
1(G) by integrating inside a predetermined algebra of polylogarithms. 

1.1.1. Three loops. There is exactly one primitive divergent graph with four loops 
[18j . namely the wheel with four spokes, or cross- hairs diagram, pictured below 
(left). Breaking it apart along each edge gives rise to exactly two topologically 
distinct bpd graphs with three loops (right). 

I ^ — T \ 2 s^~~K \ ^ — T \ 



FIGURE 2. Left: The cross-hairs diagram is the unique primitive 
divergent graph with four loops. Breaking it at the edges 1, 8 gives 
the graphs in the middle, and on the right, respectively. 



Theorem 1. Let G be one of the two bpd 3-loop graphs depicted in figure 2 (middle 
and right). Then every coefficient in the Taylor expansion of 1(G) is a rational 
linear combination of multiple zeta values. 
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Figure 3. The three primitive-divergent 5-loop diagrams. 



1.1.2. Four loop contributions. There are precisely three 5-loop diagrams, pictured 
in figure 3. The one on the left is planar and two- vertex reducible, and will be 
denoted 5 R. The one in the middle is planar, two- vertex irreducible, and will be 
denoted 5 P, and the one on the right is non-planar, and will be denoted 5 N. It 
turns out that there are exactly six topologically distinct ways to break the graph 
5 P open at one edge, and exactly two ways to break open 5 R, giving the eight 
planar topologies depicted in figure 4. 




Figure 4. The eight planar bpd topologies with four loops. 



Theorem 2. Let G be one of the eight bpd planar J^-loop graphs depicted in fig- 
ure 4. Then every coefficient in the Taylor expansion of 1(G) is a rational linear 
combination of multiple zeta values. 

For the non-planar graphs, a new phenomenon occurs, and we must introduce 
multiple zeta values at roots of unity. Let n\, . . . , n r € N, and consider the multiple 
polylogarithm function, first introduced by Goncharov [llj : 

rjM _ _ _ x k r 

(5) Lini,...,7ir( a; l) • • • > x r) = T5m Tjw ■ 

rZ-i . . . Kir 
0<fei<...<fc r 1 

It converges absolutely for \xi\ < 1 and extends to a multivalued holomorphic 
function on an open subset of C r . For m > 1, we define Z m to be the Q-algebra 
generated by the values of multiple polylogarithms at m th roots of unity: 

Lin 1 ,...,n r (afi, . • • , x r ) such that x™ = 1 for 1 < i < r , and (x r , n r ) ^ (1, 1) . 
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The condition that x r and n r are not simultaneously 1 is to ensure convergence. 
The algebra Z 1 is the algebra of multiple zeta values, and Z 2 is known as the 
algebra of alternating multiple zeta values. We will call Z m the algebra of multiple 
zeta values ramified at m th roots of unity. Note that Z a C Z h if a divides b. 




Figure 5. The three non-planar bpd graphs on four loops ob- 
tained by breaking the graph 5 N along an edge. The graph on the 
right has crossing number 1, but has been drawn with 2 crossings. 



Theorem 3. Let G be one of the three non-planar bpd graphs with 4~loops as 
depicted in figure 5. Then every coefficient in the Taylor expansion of 1(G) is a 
rational linear combination of multiple zeta values ramified at 6 th roots of unity. 

Note that the theorem gives an upper bound on the set of periods which can occur 
(Z G ), and it is an open question whether some smaller algebra Z k for k = 1,2, or 3 
suffices to compute the Taylor expansions for the non-planar graphs, or whether 
there occurs a term which genuinely involves sixth roots of unity. 

1.1.3. Five-loop contributions. Karen Yeats has computed the primitive divergent 
topologies up to seven loops [18]. There are nine at six loops, of which four are 
planar, pictured below, and five which are non-planar (figures 7,8). Although the 
five non-planar graphs all have genus 1 (they can be drawn on a torus without 
any self-crossings), it turns out that the invariant which succesfully predicts the 
outcome of our method for computing the periods at six loops is not the genus 
but the crossing number. The crossing number of a graph G is defined to be the 
minimial number of self-crossings over all planar representations of G (which are 
allowed to have curved edges). It is easy to determine the crossing number for the 
graphs we consider here, but seems to be a difficult problem in general. 

Theorem 4. Let G be a graph obtained by breaking open G at one edge, where G 
is one of the eight primitive divergent graphs at six loops pictured in figures 6 and 
7. Then the coefficients in the Taylor expansion of 1(G) are: 

(1) multiple zeta values, if G is planar (has crossing number 0). 

(2) multiple zeta values at 6 th roots of unity, if G has crossing number 1 . 

The one remaining primitive graph is the complete bipartite graph K34,, and has 
crossing number exactly 2, pictured in figure 8. Our method failed for this graph. 

1.2. Discussion. In summary, we found that for all Feynman graphs up to 6 loops, 
those with crossing number evaluate to multiple zeta values, and those with cross- 
ing number 1 give multiple polylogarithms evaluated at 6 th roots of unity. There 
is a single example of a graph (K^^) with crossing number 2, and our method fails 
in this case. The leading term in the Taylor expansion of L(G) for -£^3,4 is actually 
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Figure 6. The four planar primitive divergent graphs at six loops. 



Figure 7. The four non-planar primitive divergent graphs at six 
loops with crossing number 1. 




Figure 8. The unique non-planar primitive divergent graph at six 
loops with crossing number 2 (although six crossings are shown). 
It is the complete bipartite graph K 3> 4, and is the first graph for 
which our criterion fails. 



known to evaluate to a multiple zeta value, namely £(5,3), which, interestingly, is 
the first occurrence of an irreducible double sum in (/> 4 theory. 

One possible reason for this is that our method concerns all terms in the Taylor 
expansion of the integral 1(G), with arbitrary powers of logarithms in the numer- 
ator. In the particular case when there are no logarithmic terms in the numerator, 
the set of singularities of the integrand are slightly reduced, and it is not impossible 
that our method might work for with this restriction, although we have not 
checked this. Another possible reason for this could be because our algorithm has 
some room for improvement, in at least two ways. First, the algorithm involves 
the repeated factorization of polynomials derived from the graph polynomial. In 
our computations, we only considered factorizations that occured over the field of 
rationals Q, although it is conceivable that some of the polynomials which occur 
are irreducible over Q, but factorize over an algebraic extension of Q. 
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A second, and more promising, possibility is to extend our algorithm to deal 
with quadratic terms, and this will be explained in further detail in §4.51 It is 
possible, by extending our method in this way, that a larger class of graphs will 
become tractable. However, it is more likely than not that eventually one will find 
periods of motives which are not mixed Tate (c./. [T]), and this will pose a genuine 
obstruction to the present method. We expect that our method of polylogarithmic 
integration should also help to exhibit the first example of such a period in massless 
</>|, if and when it occurs. The idea would be to strip away from a candidate 
Feynman integral terms which evaluate to multiple zeta values, until one is left 
with a totally irreducible period integral which is vcrifiably not of mixed Tate type. 
Thus, the eventual failure of our method should at the same time exhibit the first 
non-MZV-type period, when it occurs. 

1.3. Plan of the paper. The paper is divided into two halves. The first half (§2-4) 
gives an overview of our method, and can be read linearly, as a long introduction. 
In §2 we briefly recall how to rewrite the Feynman integrals 1(G) in Schwinger- 
parametric form using graph polynomials. In §3, we outline the main idea of our 
method, and in §4, we translate this into an elementary reduction algorithm on 
graph polynomials. The main theorem [TBI in M. 41 states a sufficient condition for a 
bpd graph G to evaluate to multiple zeta values in its Taylor expansion. 

The purpose of the second half is to give a complete worked example of our 
method, in the case of the wheel with three spokes. In particular, in §7, we give 
a new proof of the well-known result that the leading term of 1(G) is 6 £(3). This 
is, to our knowledge, the first time that such a computation has appeared in print 
using a parametric representation for the Feynman integral. The purpose of sec- 
tion 5 is to provide sufficiently many details to understand the intricacies of the 
method in general, and section 6 provides worked examples of taking primitives 
and limits of poly logarithms in two variables, which are then applied directly in 
the 6 C(3) computation. It is perhaps advisable, after having read §2 — 4, to refer 
directly to the example in §7, and then read §5 and §6 bearing the example in mind. 

The proofs of our results will be written up in full detail in [9] , including a study 
of the underlying algebraic geometry, which is completely absent from this paper, 
and a generalization of the above results to some infinite families of graphs. Much 
of the background on iterated integrals, and also our method of integration with 
polylogarithms, is explained in detail in [7J [5] in the case of moduli spaces of genus 
curves. These may serve as a useful introduction to this paper, since there are many 
similarities with the integrals we study here, even though the geometry underlying 
higher loop Feynman integrals is considerably more complex. 

1.4. Acknowledgements. I owe many thanks to John Gracey and David Broad- 
hurst for explaining to me the problem of studying the three-loop contributions to 
the massless two-point function, and also to Dirk Krcimer for many helpful discus- 
sions. The list of primitive divergent graphs up to 7 loops was kindly given to me 
by Karen Yeats, and my interest in the subject owes much to the work of Bloch, 
Esnault and Kreimer [4], and the encouragement of Pierre Cartier. 

This work was motivated by the fruitful workshop on perturbative approaches 
to quantum field theory at ESI Vienna, and was continued at the MPIM during an 
IPDE fellowship, and completed at the IHES. 
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2. Parametric representations 

Let G be a Feynman graph with h loops, L internal edges, and E external legs. 
The graph polynomial of G is a homogeneous polynomial of degree h in variables 
a\, . . . , ai, indexed by the set of internal edges of G. It is defined by the formula: 

(6) u G = ^2Yl a t- 

The sum is over all spanning trees T of G, i.e., subgraphs T of G which pass through 
every vertex of G but which contain no loop. Next consider the homogeneous 
polynomial of degree h + 1 defined by: 

(7) v G = e n <* (? s ) 2 • 

s e^s 

The sum is over graphs S C G with exactly two connected components S — T± U T 2 
where both T\ and T 2 are trees, such that S is obtained by cutting a spanning tree 
5" along an edge e, and q is the momentum flowing through e in S' . 

Example 5. Let G denote the two-point two- loop diagram depicted in Fig[TJ Then 

(8) Uq = (or + a 5 )(a 2 + 04) + 0.3(011 + a 2 + a A + a 5 ) , 

Vg = (^(ot + o^Xa^ + 015) + (a2a.4&5 + Q.\a^a^ + a 1 a 2 a 5 + a 1 a 2 Q!4)) Q 2 ■ 



Figure 9. On the left are shown the eight spanning trees for G, 
corresponding to the eight terms in U G . On the right are the eight 
pairs of trees T\ U T 2 which correspond to the eight terms in Vq ■ 



To each internal edge £ of G we associate a decoration ai, which is a positive 
real number (the power to which the corresponding propagator is raised) and set 

(9) a = a£ . 

i<e<L 

The unregularised massless Feynman integral Ig{q., q, D) in dimension D, expressed 
in Schwinger coordinates ([H];[3~5], (2-36)) , is : 



a 



iTt[a+h(l-D/2)]/2 hD/2 poo poo Tf a e -l iV G /U, 

( 10 ) - TT TV / / dai - dat D% ■ 

[h r W Jo Jo Uq' 

Now let A denote a non-empty set of internal edges of G. By making the change 
of variables a.; = a', t t, for 1 < i < L, where t — a ^ integrating with respect 
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to t from to oo, and finally replacing a[ with on once again, this integral can be 
rewritten ([15]. (3.32)): 

(u) Tr hD /*T(a-hD/2) I IIW 1 fUgY-hP/2 ^ 



n,T(a,) 7 ffA c-o/a VV G V 



Here, fi L = J]f =1 (-l) l a l( iai . . . da"i . . . c?ol, and iJ A = {« t : J2ie\ a l = Tne 
fact that the integral does not depend on the choice of subset A follows from the 
fact that it is really a projective integral, but is known in the physics literature as 
the Cheng- Wu theorem. In this paper, we shall always take A to be a single edge, 
and set D = 4 - 2 e. 

2.1. Primitive divergent graphs. We say that G is broken primitive divergent 
if it has exactly 2 external legs, contains no divergent subgraphs, and satisfies: 

(12) L = 2h+1. 

In this case, let q denote the momentum entering or leaving each external leg. By 
([7]), we have Vg = Vg q 2 , where Vq is a polynomial in the Oj and does not depend 
on q. Up to gamma-factors, we can therefore rewrite (TTTj) as the integral: 

< 13 » '<--■«■< l^n(^)"'-'(f^)'. 

which converges for all Rea^ > 0. As is customary, we set 

(14) a e = l + n e e, for 1 < £ < L , 

where ng are positive integers. Before taking the Taylor expansion with respect to 
e, we first express the integrand of (|13[) as an integral of a simpler function. 

When G is broken primitive divergent, we can close up the two external legs of 
G to form a graph G with h + 1 loops, and L + 1 = 2(/i + 1) edges (see figure 1, 
right). In [4], such graphs were called primitive divergent, and one verifies that 

(15) U d = V G + a L+1 U G , 

where cxl+i is the parameter attached to the edge of G obtained by gluing the 
external legs of G together (see ^2. 21 below). The following lemma is nothing other 
than the definition of Euler's beta function. 

Lemma 6. Let < r < s, and u, v ^ 0. Then 

f°° x^ 1 , 1 T(r)T(s-r) 
/ 7 ^dx — . 

J (ux + v) s u r v s r r(s) 

On comparison with (jlip . we set r = (h + l)-j — a, and s = -j. Substituting 
([T2"]) and (fT4"]) into the expression for r leads us to define: 

L 

(16) a L+1 = 1- (ft + l + ^ni)e . 

»=i 

Therefore, using the previous lemma, we can rewrite the integral (|13p as the product 
of certain explicit gamma factors with the following integral: 

r n i+1 a a4 " 1 

(17) I( ai ,...,a L+1 ,e)= lU %_i n I+ i, 
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where Ol+i = H,LdaL+i, and H\ is the hyperplane a\ = 1 in [0, oo) L+1 . The 
Taylor expansion of the original Feynman integral (llip can be retrieved from the 
following well-known formula for the Taylor expansion of the gamma function: 



oo 

T(s + 1) = exp(- 7S ) exp ( ^(-l) fe C(fc)s' 



k=2 

The coefficients in the Taylor expansion of (|17[) are given by: 

(18) / ^P(ai, log(oi), log(CTg)) , 

where m is an integer, and P is a polynomial in oti,log(oii) for 1 < i < L + 1 and 
logt/g with rational coefficients. This is a period integral in the sense of [T5] , 

In the first half of this paper, we will state a criterion on G for (fTSf to evaluate 
to multiple zeta values, and in the second half, we will outline an algorithm for 
computing these values. 

2.2. Contraction and deletion of edges. For any graph G, its graph polynomial 
Uq is linear with respect to each variable at. We can write 

(19) U G = U G//{i} + f/ G \ {l} a l , 

where G//{i} is the graph obtained from G by contracting the edge labelled i, and 
G\{i} is the graph obtained from G by deleting the edge labelled i. Using (fTS)) . 
this implies that 

(20) U G = U GX{L+1} and V G = U d//{L+1} . 
From (fT9f and (j20|) . one deduces a similar formula for Vg, namely: 

(21) V G = V G//{i} + Vcx^a, . 

Notation 7. When the graph G is implicit, we will often write U, V instead of 
Ug,Vg- In order to lighten the notation, we adopt the useful convention from [3], 
which consists in writing, for any polynomial ^> which is linear in the variable c^: 



d 

<I> W = - — * , and *j = * 

OCti 



a i= 



We therefore have Ut — Uq/iux, U^ l > = U G \{i}, and (fTl?)) and (|2T|) become 

17 = Z7i + U il) a t and 7 = + . 
Since the operations of contracting and deleting distinct edges commute, we can 
write U12 = U G //{i,2}, U% = U G //{2}\{1}, U {12) = C/ G \{i,2}, and so on, where 
indices in the superscript (subscript) correspond to deleted (resp. contracted) edges. 
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3. Symbolic integration using polylogarithms 

Let G be a broken primitive divergent graph, and let G be the graph obtained 
by closing its external legs. To illustrate our integration method, let us begin to 
compute the integral (fl3|) in the case where all decorations ai — 1, and e = 0. 
Therefore consider the convergent integral 

(22) I S = f %±i , (Step 0) 

which was studied in [4]. We will assume that A is a single index, say L. In this 
case we can write £Il = da\da 2 ■ ■ ■ rfa^-i and = ^1^1+1- The domain of 

integration is simply H\ = {o>l = 1, < a* < 00 , i ^ L}. Using equation (fT5"l) 
to replace t/g with Vg + c*l+\Ug, we can perform one integration with respect to 
oil+i from to 00, to obtain: 

da L+1 n f tt L 



(23) Ig = (v - , 2 V L = / -±- (Step 1) 

{V G + a L+ iU G y J Hx UgVg 



From here on, we write U, V instead of Ug, Vg and use the notations of W2.'2\ Now 
choose any variable in the integrand, say a±, and write Iq in terms of a\\ 



Ic^iiUt + UWaMVt + VWa!) 
By decomposing into partial fractions with respect to ai, this gives 

i G = 



Hx U^Vi - UxVW \Ux + UWat Vi + V^ai 

Now we can peform an integration with respect to u\ from to 00 at the expense 
of introducing a new function, namely the logarithm. 

(24) Ig = I ^-^-^V^^ da2 . . . daL _ x (step 2) . 
v ; G J aL=1 I7(i) Vi _ UlV (i) L 

The logarithm log / should be regarded as a symbol which satisfies the formal rule 
d\ogf = f^ x df, since changing the constant of integration does not affect the 
integrand of (|24| . At this stage, something remarkable occurs. The denominator 
factorizes as the square of a polynomial D which is linear in each variable 0:2, • • • , oil: 

(25) D 2 = U {1) V! - £/iU (1) 

As observed in ,4], this phenomenon is quite general and follows from a result due 
to Dodgson on determinants of matrices. One can also give a formula for D in 
terms of trees in G, but this will not be required here. 

We can repeat this argument. Choose a variable, say ct2, and write 

D = D 2 a 2 + 7J (2) . 

By decomposing the integrand of (|24|) into partial fractions with respect to the 
variable a 2 , we can integrate with respect to the variable a 2 - To simplify the 
notation, we define 

, oa , r 1 , p(log(p) +log(s) -log(g) -log(r)) log(g) 

26 {p,qr,s\ = 1 . 

r(ps — qr) rs 
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One verifies that 

Corollary 8. The integral Iq is equal to the L — 2 dimensional integral: 

Ja L =l V 

+{Vi2,vl 2) \D2,D^}^da 3 ...da L ^ (Step 3) . 

As long as there exists a variable with respect to which all polynomials which 
occur in the integrand are linear, this can be repeated. 

At the next stage of the integration process, one has to introduce the diloga- 
which is formally defined by the differential equation 

(27) dU 2 (f) = - l ^fl d f. 

In this manner, we obtain a conditional algorithm for computing integrals such 
as (I18p . which can be approximately formalised as follows: 

(1) Choose a variable in which all terms of the integrand are linear. 

(2) Formally take a primitive of the integrand with respect to this variable. 

(3) Evaluate this primitive at and oo, and repeat. 

The algorithm fails exactly when in (1), we can no longer find a variable with 
respect to which all terms in the integrand are linear. This will be formulated pre- 
cisely in the following section. It can be checked in advance whether the algorithm 
will terminate, as this depends only on the topology of the graph G. When the 
algorithm does terminate, any integral (|13[) with arbitrary decorations can always 
be computed in a finite number of steps in terms of a fixed differential algebra of 
polylogarithms determined in advance by the topology of G. 

For example, in the above, we can keep track of the singularities of the integrand 
(i.e., terms occuring in the denominator, or arguments of the logarithm functions) 
at each step. At the first step, we represent the singularities as the set: 

{U,V}. 

After the second integration, the singularities are given by: 

{uV>,U u vV\VuD} . 

At the third step (corollary [5]) , we have in the same manner 

{U^ , UP , , U* , , vf\ ,V 12 ,DW,D 2 }, 

along with the set of irreducible factors of the four denominator terms: 

D^U^-D 2 U^ , D^U 12 -D 2 u[ 2) , D^V 2 (1) -D 2 V^ , D^V 12 -D 2 vf ] . 

We are therefore led to consider an algorithm for the reduction of sets of polynomials 
with respect to their variables. This will ultimately lead to a criterion for the 
computability of broken primitive divergent Feynman graphs (theorem 1 18[) . 



As we shall see later, it is more convenient to use the function satisfying the differential 
equation dL(f) = / _1 log(/ + 1) df, which is a close relative of the dilogarithm. 
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4. Reduction of polynomials 

Keeping track of the polynomials, or singular loci, which can occur in this inte- 
gration process gives rise to a reduction algorithm which can be used to check the 
outcome of the integration process without actually doing it. 

4.1. The simple reduction algorithm. Let 5 — {/i, . . . , /a}, where /i, . . . , /jv 
are polynomials in the variables a%, . . . , a m , with rational coefficients. 

(1) Suppose that there exists an index 1 < r < m with respect to which every 
polynomial /i, . . . , /jv is linear in the variable a r . Then we can write: 

fi = gia r + hi for 1 < i < N , 

where gi = dfi/da r , and hi = fi\ ar =o- Define a new set of polynomials: 

S{r) = {(3i)l<i<AT > ( h i)l<i< N ' ( h i9l -& h i)l<i<j<N} ' 

(2) Let Sm be the set of irreducible factors of polynomials in Sm. 

The polynomials now occurring in 5( r ) are functions of one fewer variables, namely 
ai, ... , a r _i, a,+i, . . . , a m . This process can be repeated. If at each stage there 
exists a variable in which all polynomials are linear, we can proceed to the next 
stage. This gives a sequence of variables and a sequence of sets 

(28) ^(n) : S(r 1 ,r 2 ) > • • ■ ! ^(ri,...,r») 1 

If there exists a sequence (7*1, . . . ,r m ) such that every variable is eventually elimi- 
nated, then we say that the reduction terminates. When this happens, we say that 
the set 5 is simply reducible. 

Remark 9. We can remove any constants, and any monomials of the form on which 
occur as elements in Sr r ), as this does not affect the outcome of the algorithm. 

Definition 10. Let G be a broken primitive divergent graph, and let 

So = {U d } . 

We say that G is simply reducible if Sq is simply reducible. 

Observe that {Uq}(l+i) = {Ug,Vg}, where ctL+i is the edge variable of G 
obtained by closing the external legs of G, by equation (115|) . 

Example 11. Consider the 2-loop 2-point graph G depicted in fig. 1. We write 
S = Sq — {Uq}- After reducing with respect to olq, we have 5(6) = {U, V}, where 
U = U G , V = Vg are given by ©: 

5(6) = {aia2 + ol\ol^ + 015012 + 0405 + a 3 a\ + a^,a 2 + 0304 + 03015 , 

030404, + Q!3Q!iq;5 + a30!2a!4 + a3a2a5 + a2a4C*5 + aia4a5 + aio;2Q!5 + aia^a^} ■ 
Since both polynomials are linear in ai, we can reduce with respect to ot\ to obtain 

5(6,!) -{C/ (1) ,C/i, v^ 1 \v l M l) V4-u l v^} . 

But U^Vi-UiV^ factorizes, by the Dod gson identity, and so at the second stage 
we have 5( 6 ,i) = {C/ (1) , U u Vi, D}, where D 2 = (U^Vi - l^W): 

5(6,1) = {a2«5 + 0405 + 0:20:3 + 0304 + 0305 , Q-2 + O3 + a 4 ) 0304 + 0305 + 0405 , 
030:4 + 0305 + 0405 + 0205 + 0204 , 0304 + 0405 + 0205 + 0305 } 
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Since each polynomial is linear, we can reduce with respect to the variable ct2'- 
S(6,i,2) = {«4 + as , a 3 a 4 + a 3 a 5 + a 4 a 5 , a 3 + a 4 , a 3 + a 5 , a 3 - a 4 } . 
Next, reducing with respect to the variable 0:5, gives 

5(6,1,2,5) = ( a 3 + a 4 , a 3 - a 4 } , 
and finally, reducing with respect to a 3 gives 

5(6,1,2,5,3) = • 

Here, and from now on, we will adopt the convention (remark [9} that we remove 
all constant terms and variables on from each stage of the reduction. 

In <j3l the domain of integration for Iq was a hyperplane a\ = 1 for some index A. 
If Sq is simply reducible with respect to some order (a ri , . . . , a Tn ) of the variables, 
we can choose A to be the index of one of the two final variables In the 

previous example, we can take A = 4, and set a 4 = 1, S' = SG\ a =1 - This gives 

5(6,i,2) = {"5 + 1 , a 3 + a 5 + a 3 a 5 , a 3 + 1 , a 3 + a 5 , a 3 - 1} , 

and 

5(6,1,2,5) = {"3 + 1 > "3 - 1} • 

The reduction algorithm reflects the set of singularities which occur at each stage 
of the integration process. Roughly speaking, the integrand at the fc th stage, after 
integrating with respect to a ri , . . . , a rk , will have singularities along the zero locus 
of polynomials in the sets 5( ril ... lT - m ), for m > k, and along the axes on = 0. 

Remark 12. The algorithm described above is similar to an algorithm defined by 
Stembridge to study the zeros of graph polynomials over finite fields. A similar 
argument was also used in [1] to obtain a Tate filtration on graph hypersurfaces. 

In the previous example, we see that at the penultimate stage, we expect to have 
polylogarithms in a 3l with singularities along a 3 = 0, a 3 = — 1, a 3 = 1. This would 
imply that at the final stage, the integral Iq is a rational linear combination of 
alternating multiple sums iJ 2 , i.e., periods of P 1 \{0, — 1, 1, 00}. This is not good 
enough, since it is known that the integral Iq in fact gives multiple zeta values. To 
rectify this problem, we introduce the Fubini reduction algorithm. 

4.2. The Fubini reduction algorithm. In the simple reduction algorithm, the 
sets S( ri rt,) which control the singularities of the integration process, depend in 
an essential way on the order of the variables n, . . . , 7> which was chosen. However, 
it is an obvious consequence of Fubini's theorem that the final integral does not 
depend on the particular order of integration. More precisely, if F is the integrand 
obtained at the fc th stage, we clearly have 

p 00 />oo p 00 /*oo 

/ da ri / da rj F = / da r . / da Ti F . 
Jo Jo 3 Jo 3 Jo 

The left-hand integral has singularities contained in the zero locus of polynomials in 
5(ri,...,r fc ,r 4 ,r-)5 the right-hand integral has singularities contained in the zero locus 
of polynomials in St ri rk>rj ri y It follows that both have singularities contained 
in the zero locus of polynomials in the intersection 5( I . 1) ... jrfciT . ijr .) H S'(r 1 ,...,r fc ,r ;j ,r i )- 
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In general, we define sets recursively as follows: 

(29) <5Vi,r 2 ] = S(n,ra) n ^(ra,n) > 

S[ ri ,r2,...,r k ] = ^lri,...,fl,...,rk](n) > ^ > 3 , 

l<i<fc 

where in the second line one applies a reduction to the set S , r ri)- . -i fj ) . - . iT .j with respect 
to the variable a Ti (steps (1) and (2) in the previous section). It may happen that 
it is not possible to apply the reduction to S\ ri ,...,r k ] if it contains a polynomial 
which is non-linear in a Ti . In that case, S , [n....,f~ i ,...,r fc ](r i ) is undefined and we omit it 
from the intersection in the second line of |29|) . It is also possible that this happens 
for all 1 < i < k, in which case S[r\, . . . , rk] is undefined. 

This is the Fubini reduction algorithm, and it gives rise to a sequence of sets 

(30) ^(ri) > S[ri,r 2 ] > • • ■ ! &[ri ,...,r n ] ■ 

Definition 13. We say that S is Fubini reducible if there is a sequence (ri, . . . , r m ) 
such that every variable is eventually eliminated, and such that every polynomial 
in S\ rii ... t r k ] is linear in av fc+1 , for all 1 < k < m. 

If G is a broken primitive divergent graph, we say that G is Fubini reducible 
if Sq = {Uq} is Fubini reducible. As in remark [9l we can remove constants and 
terms of the form on without affecting the algorithm. 

Remark 14. By the definition (I29p . the set S\ ri>r2t ,,, irk ] is contained in the set 
5'[r T(1) ,r^ (2) ,...,r„ (fc) ] for any permutation n of 1, 2, . . . , k. Thus if S is simply reducible, 
then the sequence ([30|) is contained in the sequence ([28)) . In general, the inclusion 

"5 , [ri,r 2 ,...,r fc ] £ ("j £(7r(ri),...,ir(rfc)) j 

7ree fc 

is strict, where the intersection is taken over all permutations of ri, . . . , r&. In 
practice, the cascade of polynomials one obtains using the Fubini reduction is con- 
siderably smaller than with the simple reduction algorithm. 

Example 15. We compute the Fubini reduction algorithm for the two-loop two 
point function, as in example [TT] Let S — {Uq}, where Uq = ocq Ug + Vg, and 
Ug, Vg are given by ([8]). One checks that Spj.eirs) is given by: 

{«3+a4 ,a3+a4+ai , axa^+aia^+azct^ ,ai+a3 ^a^a^+o^+aiai+aia^ , aia^+aia4+a^} 
and 5'[5,6](2) by 

{03+04 ,ct3+a4+ai , aia4+Q;ia3+Q!3a4 ,ai+a^ , aia4+o;3+a3a4 , 2aia3+aick4+a!3+a3a4} 
and 5[ 2 , 5 ](6) by 

{a 3 + «4 , a3 + a4 + ai , ol\ol± + aia 3 + a^a^ , ot\ + 03 , a± + 0:4, a\ — a^} . 
Taking the intersection of all three sets gives 

^[2,5,6] = i a 3 + a ± i a 3 + a 4 + a l , "1^4 + aiC*3 + a 3 0!4 , Oil + a 3 } . 

By performing an ordinary reduction with respect to the variable a%, one obtains: 
S[2,5,6](i) = {«3 + «4, a z a± + al + al} . 

By computing SV 6 ,1,21(5)) £[6,1,51(2)! £1,2,51(6) and intersecting all four (actually it 
suffices to intersect with 6(6,1,2,5) from example I 111 in this case), one verifies that 

£[6,1,2,5] = {"3 + Oii} . 
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As before, we can take our domain of integration to be the hyperplane 04 = 1. At 
the final stage of the integration process, the Fubini reduction algorithm will predict 
(see theorem [T5| that the integrand has singularities in 013 = 0, a 3 = 1. Therefore 
we expect to obtain multiple polylogarithms with singularities in {0, 1, 00} at the 
penultimate stage, and hence multiple zeta values as the final answer. In this way, 
the problem raised at the end of the previous section has apparently been overcome. 

However, constants can appear at every stage during the integration process (see 
§7), and one needs to verify that we only obtain multiple zeta values every time. 
This leads to a further ramification condition to be verified, for each set S\ rij .„ trk ]. 

4.3. The ramification condition. Let S be Fubini reducible for some order 
(7*1, . . . ,r m ) of the variables. After setting the final variable a Tm = 1 in the se- 
quence ([30]) . we obtain a new sequence of sets (the reductions of S' = S\ a 

(31) ^Vi) ' ^[ri,r a ] ' ■■■ , »S'[r 1 ,...,r m _i] > 

where every polynomial in SL r ^ is linear in the variable ct rk+1 . Therefore, if we 
write S' [ru rk] = . . . , fM k }, then we have 

fi = at a rk+l + h , for 1 < i < M k , 

where a,, bi are polynomials in a rfc+2 , . . . , a Tm . We define T, ak to be 

(32) S Qfc = {-— such that a, ^ 0} . 

a l 

The set S Qfc clearly depends on the ordering of the variables (ri, . . . , r m ). 
Definition 16. We say that S Qfc is unramified if: 

(33) lim ( lim (...( lim • • •) C {0, -1, 00} . 

We say that the sequence of sets (|3"0)) is unramified if the corresponding S af! are 
unramified for all 1 < k < m — 1 . Finally, we say that S is unramified if it is Fubini 
reducible, and if there exists a sequence (n, . . . ,r m ) such that (|3"Tj) is unramified. 

Example 17. One can check that the sequence of sets given in example fTS] is 
unramified. We check the last two terms only. Setting = 1, we have: 

^[6,2,5] = {1 + "3 , l + ai+a 3 , ax + a 3 + «ia 3 , a.\ + a 3 } , S[ 612 5] = {1 + a 3 } . 

Thus 

S ai = {-(ata + 1) , -— , -a 3 } , and E Q3 = {-1} . 

1 + a 3 

The ramification condition (|33p is satisfied in both cases, since Eq, 3 C {0, —1, 00}, 
and lim Q3 ^o S ai = {0, —1}. 

4.4. The main theorem. The main theorem gives a simple criterion for a master 
Feynman integral to evaluate to multiple zeta values. 

Theorem 18. Let G be a broken primitive divergent Feynman diagram and let 
Sg = {Uq}- If Sg is Fubini reducible and unramified, then the coefficients in 
the Taylor expansion of the integral 1(G) ( |_?7[ ) are rational linear combinations of 
multiple zeta values. 
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We also require a variant to allow for multiple zeta values ramified at roots of 
unity. Let S be Fubini reducible, as above. We say that S is ramified at p th roots 
of unity if there exists a sequence (r±, . . . , r m ) such that the corresponding sets S afc 
satisfy 

S Qfe C {0,oo}U{-cj:w p = 1} . 
If p = 1, then this coincides with the definition [TcjI above. 

Theorem 19. Let G be a broken primitive divergent Feynman diagram and let 
Sq = {Uq}- If Sq is Fubini reducible and ramified at p th roots of unity, then 
the coefficients in the Taylor expansion of the integral 1(G) are rational linear 
combinations of multiple zeta values ramified at p th roots of unity, i.e., in Z v . 

The results stated in the introduction are the result of applying these two the- 
orems to the set of all primitive divergent graphs up to 5 loops. We computed 
the Fubini reduction algorithm for the cross-hairs diagram (for which the graph 
polynomial has 45 terms), the graphs 5 R, 5 P and 5 7V (for which it has 128,130 
and 135 terms, respectively), and all eight primitive divergent graphs with 6 loops 
and crossing number and 1, as depicted in figures 6 and 7 (for which the graph 
polynomials have around 400 terms). We found that all planar bpd graphs up to 
5 loops are Fubini reducible and unramified, and all non-planar bpd graphs with 
crossing number 1 are Fubini reducible but ramified at 6 th roots of unity. The 
bipartite graph K^^ has crossing number 2 (figure 8), and is not Fubini reducible. 
These calculations were done using maple. 

The outcome of the Fubini reduction algorithm may depend on the particular 
choice of two final variables. Typically, the ramification is better if one chooses the 
two final variables to correspond to edges meeting at a four-valent vertex. In each 
example, we computed the algorithm for all possible pairs of final variables, and 
chose a pair giving the least ramification. In an ideal world, one could simply take 
the intersection of the period rings obtained for every such choice, and the algorithm 
would in that case not depend on any choices. In general, however, taking the 
intersection of rings of periods requires some powerful diophantine results, which 
are at present totally out of reach. 

4.5. Extensions. There are a number of ways in which the reduction algorithm 
might be improved. First of all, in our computer calculations, we factorized our 
polynomials over Q, rather than over the algebraic closure Q, which may or may 
not have made a significant difference. 

More interestingly, the Fubini reduction method stops as soon as it finds that, 
for every variable on, there is a polynomial which is quadratic in that variable. It is 
possible, by introducing new variables which are the square roots of discriminants 
(by passing to a ramified cover), one could generalize the method to deal with some 
of the quadratic terms. The basic idea is that, in the case of a plane curve for 
example, a polynomial need not necessarily be linear to define a curve of genus 0. 
A conic, of degree 2, also defines a curve of genus 0, and our method should also 
extend to this case. A genuine obstruction should occur at the degree 3 level, since 
one expects to find elliptic components in this case. In any case, by extending our 
method to deal with quadratic terms, a larger class of graphs may then become 
reducible, and amenable to computation. 

Finally, if one is interested only in the leading term of the Taylor expansion of 
the graph 1(G), then one can simplify the reduction algorithm. Recall that, for the 
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general integral (with arbitrary logarithms in the numerator), we obtained at the 
second step, the set of singularities: 

S [L+ltl] ={uW,Ui,vW,V u D} . 

Notice that in (f!M)) . where the numerator has no logarithmic terms, not all pairs of 
singularities can occur. Thus one should perform a Fubini reduction separately on 
the four sets 

{UU,D} , {U U D} , {V^\D} , {V U D} , 

which may improve the ramification. There may also exist graphs which are not 
Fubini reducible, but for which each of the four sets above is Fubini reducible. This 
would prove that the first term only in the Taylor expansion is a multiple zeta value. 

In the remainder of the paper we explain how to compute the coefficients in the 
Taylor expansion when the conditions of theorems [T8l and fT9l hold. 
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5. HYPERLOGARITHMS, POLYLOGARITHMS, AND PRIMITIVES 

In this section, we outline the function theory underlying the iterated integration 
procedure. For a more detailed introduction, see [51 El IB] • 

5.1. Hyper logarithms. Let £ = {(To, &i, ■ ■ ■ > cat}, where «7j are distinct points of 
C. We will always assume that <j$ — 0. For now we will consider the points cr; to 
be stationary, but later we will allow them to be variables in their own right. 

Let A = {ao, ai, . . . , a at} denote an alphabet on N+ 1 letters, where each symbol 
3i corresponds to the point <7j. Let Q(A) denote the vector space generated by all 
words w in the alphabet A, along with the empty word which we denote by e. To 
each such word w, we associate a hyperlogarithm function: 

L w (z) : C\S — ► C 

which is multivalued, i.e., it is a meromorphic function on the universal covering 
space of C\S. Let log(z) denote the principal branch of the logarithm. 

Definition 20. Let A x denote the set of all words w in A, including e. The family 
of functions L w {z) is uniquely determined by the following properties: 

(1) L e (z) = 1, and L a n(z) = ± log"(z) , for all n > 1. 

(2) For all words w £ A x , and all < i < N, 

d 1 

^-L 3iW (z) = L w (z) , for z G C\S . 

oz z — Oi 

(3) For all words w G A x not of the form w = ag , 

lim L w (z) = . 

The weight of the function L w (z) is defined to be the number of letters which 
occur in w. The functions L w (z) are defined inductively by the weight: if L w (z) 
has already been defined, then L aiW (z) is uniquely determined by the differential 
equation (2) since the constant of integration is given by (3). 

It follows from the definitions that 

L 3t (z) = log(z - Gi) - logO;) . 

Later, we will only consider linear combinations of functions L w (z) which are single- 
valued on the real interval (0, oo). Any such function is uniquely defined on (0, oo) 
after having fixed the branch of the logarithm function L ao (z) — log(z). 

5.1.1. The shuffle product. The shuffle product, denoted m , is a commutative mul- 
tiplication law on Q(A) defined inductively by the formulae: 

wme = eniw — w for all w G A x , 

a;u>i majW2 = 3i(wi ma^u^) + aj(aiWi niW2) for all w\, W2 G A x , a,;, Sj E A . 
We extend the definition of the functions L w (z) to Q{A) by linearity: 

m m 

L w {z) = qiL Wi (z) where w = q i w l , q { G Q , u>i G A x . 

i=l i=l 

The following lemma is well-known. 

Lemma 21. The functions L w (z) satisfy the shuffle relations: 

L Wl (z)L W2 (z) = L WllJIW2 (z) , for all wi,w 2 G Q{A) , and z G C\E . 
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Definition 22. Now let us define 

1 1 



0E=' 



Z — tjQ Z — (TJV- 

and let L(E) be the 0£-module generated by the functions L w (z), for w G A x . 

The shuffle product makes into a commutative algebra. It is a differential 

algebra for the operator d/dz, and is graded by the weight. 

Theorem 23. The functions {L w (z)} we Ax , are linearly independent over C<g>Q0£- 

It follows from the theorem that is a polynomial ring, and a convenient 

polynomial basis is given by the functions L w (z), where w are Lyndon words [14] . 

In order to take primitives in we have to enlarge the ring of coefficients 

slightly. Therefore we must consider 



at, " , < i < j < N 



and let L + (T,) = 0^ ®o-e be Oj-module spanned by L w (z), for w G A x . 

Theorem 24. Every element f in L(E) of weight n has a primitive in L + (E) of 
weight at most n + 1, i.e., an element F G L + (S) such that 

dF/dz = f . 

The primitive of any generator f(z)L w (z) G where f(z) G Oe, can be 

found explicitly by decomposing f(z) into partial fractions in O^. The fact that 
one must enlarge Os to 0\ is clear from the identity: 

1 1/1 1 



[z - CTj)(z — (Tj) <Ji - o 

One is thereby reduced to the case of finding a primitive of functions of the form 

(z — ai) n L w (z) , where n G Z . 

In the case n = —1, a primitive is given by L 3iW (z) by definition (2) above. In 
all other cases, integration by parts enables one to reduce to the case of finding 
a primitive of a function of lower weight. Thus a primitive of f{z)L w {z) can be 
computed algorithmically in at most n steps, where n is the weight of L w (z). 

5.1.2. Logarithmic regularization at infinity. Having fixed a branch of the logarithm 
L ao (z) we can define the regularization of hyperlogarithms at infinity. 

Proposition 25. Every function f{z) G L(E) can be uniquely written in the form 

m 

f( Z ) = Y J m^g\z) , 

where fi(z) is holomorphic at z — oo, for < i < m. 

We can therefore define the regularized value of / at infinity to be: 

Reg 2=oc /(z) = /o(oo) . 

Clearly, Reg z=oc L ao (z) = 0, and for all 1 < i < N, Reg z=oc L ai (z) = -log(tJi). 
Note that the regularization operator Reg z=OQ respects multiplication. 



THE MASSLESS HIGHER-LOOP TWO-POINT FUNCTION 



21 



There is an analogous notion of regularization at z = 0. Note that by the 
definition of the functions L w (z) (properties (1) and (3)), we always have 

(34) Reg z=0 L w (z) = , for all w £ A x , w ^ e . 

Definition 26. Suppose that f(z) £ L(E) is holomorphic on the real interval 
(0, oo), and that f(z)dz has at most logarithmic singularities at z = 0, oo. We 
define the regularized integral of f(z) dz along [0, oo] to be: 

y* oo 

/ f(z) dz = Reg z=OQ F{z) - Reg z=Q F{z) , 
Jo 

where F (z) £ L(Y,) is a primitive of f{z). The integral converges. In practice, we 
write F(z) in terms of the basis of functions L w (z), and choose the constant term 
to be zero. It will follow that Reg z=0 F(z) vanishes, and the integral is simply given 
by Reg z=OQ F(z) in this case. 

It is clear that the regularized integral is additive. The point is that, in order to 
compute an integral which is convergent, one is allowed to break it into a sum of 
logarithmically divergent pieces, compute the regularized integral of each, and add 
the answers together. This is illustrated in ^7.21 

5.1.3. Multiple Zeta Values. Consider the case N = 1, and ctq = 0, a% = — 1. After 
a change of variables z \— » —z, we retrieve the classical situation on P 1 \{0, 1, oo}. 
The functions L w (—z) are therefore multiple polylogarithms in one variable [7], and 
we can compute the regularized values at infinity in terms of multiple zeta values 
using standard methods. To do this, let X = {xo,Xi} denote an alphabet with two 
letters, and let Z(xo,Xi) be Drinfeld's associator: 

Z(x ,xi) = ^2 Cm{w) w , 

ujexx 

where the numbers £hi (w) £ R satisfy: 

Cm (xo) = Cm (xi) = , and Cm (xq^Xi . . . Xq^Xi) = £(m, . . . , n k ) , 

for all integers m, . . . , njt > 1 such that n k > 2. Extending Cm (w) by linearity to 
Q(X), we have the shuffle relations 

Cm (wi)Crn {w 2 ) = Cm (»i mw 2 ) for all wi , w 2 £ X x . 

The numbers Cm (w) are uniquely determined by these properties. Using the fact 
that C(l, 2) = C(3), one can verify that, up to weight three: 

Z(x ,xi) = 1+C(2) (x Xi-XiX )+C(3) (xoXi+x Xi+XiXo+XiXo-2xoXix -2x 1 xoXi)+. . . 
Lemma 27. Let S = {0, — 1, oo}, and let w £ X x . Set 

(35) CooH = Reg z=OQ L w (z) . 

Then the generating series of regularized values at infinity is given by: 
^2 C<x>(w)w = Z~ 1 (x 1 -xo,-xo) . 

A simple computation shows that: 

(36) Z- X (xi -x , -xo) = 1 + C(2) (xix - xqXi) + C(3) (x x? + x?x - 2 xix xi) + . . . 
which allows us to read off the values of Coo{w) for w up to weight 3. 
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5.2. An algebra of polylogarithms for reducible graphs. The integration 
process takes place in an algebra of polylogarithm functions in several variables. 
A basis of such functions will be products of hyperlogarithms L W1 (z±) . . . L Wk (z k ), 
each of which is viewed as a function of a single variable Zi, where z% is some 
Feynman parameter. Now, however, the singularities oq, . . . , on will themselves be 
rational functions in the remaining Feynman parameters ax, . . . , akO 

Let S — {fx, . . . , /m} denote a Fubini reducible set of polynomials in the vari- 
ables ai, ... ,ajsf- We can assume that the algorithm of H4.'2\ terminates when we 
reduce S with respect to the variables ax, . . . , Ozv in that order. We define a nested 
sequence of rings Rn C Rn-1 C . . . C Rx recursively, as follows. Set Rn+i = Q 
and define: 

(37) R k = R k+1 [a k ,— ,~ :/eS [ai , aj at| ] , for 2 < k < N . 

OLk J 

The fibration algorithm ensures that every / 6 ST ai)a2 . iafc i is linear in a k - Thus, 
if {fx, . . . , f mh } is the set of elements in Sj aija2) ,„ )afc ] whose coefficient of a k is 
non-zero, we can write: 

fi = at a k + bi , where a,, bi are invertible in Rk+x , 

for all 1 < i < rrik- As in i )4.2i we set 

(38) E fc = {ax, ■■• ,cr mfc } , where a t = -— , 

ai 

and with these notations, we can write: 

(39) R k = Rk+x \a k , , en G Sfc] 

ak - Oi 

We refer the reader to §7.11 for a worked example in the case of the wheel with 3 
spokes diagram. 

Definition 28. Let L(R]y) — R^, and define inductively 

(40) L(R k ) = i(Sfe) ® L{R k +x) iovl<k<N . 

Every element of L(Rk) can be represented as a sum of terms of the form: 

(41) 4> = f(a k , ■ ■ ■ ,a N ) L Wl (ak)L W2 (a k+ x) ■ ■ ■ L WN k+1 (a N ) , 

where / € Rk is a rational function of a k , ■ ■ ■ , (X-n, an d L Wi (&k+i-i) G L(J^k-i+i) is 
a hyperlogarithm in a k +i-i, whose set of singularities Y.k-i+x are rational functions 
of the higher variables ctk+i , ■ ■ ■ , «at ■ 

Definition 29. The weight of the element <fi is |wi | + . . .+ \ wn\, where \w\ denotes 
the number of letters in a word w. 

Every element 4> S L(Rk) is a multivalued function on a certain open subset of 
C N ~ k . When we integrate, we will only consider elements (f> which are holomorphic, 
and hence single-valued, on the real hypercube (0, oo) N ~ k C C N ~ k . 



We use the word hyperlogarithm to denote a function L w (z), considered as a function of the 
single variable z, with constant singularities cr^ £ C. In contrast, we use the word polylogarithm 
to denote L w (2) also, but this time viewed as a function of several variables 2 and , . . . , ojy 
where some of the singularities ai depend on the a k , . ■ . ,aty 
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In this manner, we have defined an explicit algebra of poly logarithm functions as- 
sociated to a reducible Feynman graph G. Its elements can be represented symbol- 
ically as a linear combination of products of words. We next show how to compute 
any integral of the type ([4]) associated to G by working inside this algebra. 

5.3. Existence of primitives. The integration process requires finding primitives 
at each stage. As explained earlier, we need to enlarge the coefficients of the algebra 
L(Rk) slightly in order to do this. In the above notation, we set 

Rfc = Rk &i, , for <Ti,<Tj G Sfe distinct . 

0~{ 0~j 

The following theorem establishes the existence of primitives. 



and Reg ak=Q F = 



Theorem 30. Let f G L(Rk) denote a function of weight n in the variables 
ctk, ■ ■ ■ ,ctisf. Then there exists a unique function F S R^ ®R k L(Rk) of weight 
at most n + 1, which is a primitive of f with respect to ctk, and is regularized at 0: 

OF 

dak 

The construction of the primitive follows immediately from theorem [Ml We can 
assume / is a generator of L(Sfc) of the form / = fig, where f\ £ L(Y.k), and 
g G L{Rk+\) is a function of a k +i, ■ ■ ■ , a^ only. Theorem 1241 provides a primitive 
F\ G R~l ®R k L(Rk) for f\. The required primitive of / is therefore F — F\g. 

5.4. Restricted Regularization. In the notations of theorem [301 we must next 
consider the regularized limit of such a primitive F at a^ = oo, 



r-OG 

/ f da k = Rcg ak=oc F . 
Jo 



This is a function of a k +i, ■ ■ • , a at, but will not lie in C(afe_|_i, . . . , ajv) ®q L(Rk+\), 
since it will have extra singularities corresponding to the loci oi = o~j , where &i , o~j G 
S/j. The main technical point is that these extra singularities will cancel during the 
iterated integration process. Therefore, by mapping superfluous singular terms to 
zero, we can compute the integrals by using a more economical space of functions. 
In this way, we will define a restricted regularization map: 

(42) RReg afc=oc : L(R k ) — > C{a k+1 , . . . , a N ) <8> Q L{R k+x ) , 

which consists of first taking the regularization at infinity as defined in H5.1.2I and 
then projecting onto the algebra C(afe-|_i, . . . , aw) ®q L(Rk+i). 

The projection map is easily illustrated for functions of a single variable, i.e., in 
the case of hyperlogarithms. Consider two sets of distinct points in C: 

S = {cro,o'i,...,o'„}Gi; / = {do, ■ • ■ , cr n ,a n+ i, . . . , a m } , 

and let L(T,') denote the corresponding hyperlogarithm algebras, indexed by 

the set of words in the alphabets A = {a , . . . , a„} C A' = {a }, respec- 

tively. There is a projection map 

(43) 7T S :£(£') — ► Oe<®q£(£) 

L w {z) i > L 7rA{w) (z) , 

where ka{w) = w if w G A x , but ka{w) = if w contains a letter in A'\A. The 
map 7T£ is a homomorphism for the shuffle product. Thus, for any set £' containing 
S, we can in this way project out any singularities in S' which are not in S. 
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The idea is to apply this argument to Reg afc=DO F. If one considers afc+2, ■ ■ ■ , ctN 
to be constant, then one can verify that Heg ak=00 F is a hyper logarithm function 
in the variable ak+i- We can therefore write it in the form: 

Reg ak=oc F = ^qi L Wi (a k +i) , where L Wi (a k+1 ) G £(£') , 

i 

where L Wi (ak+i) has singularities in some set £' which may be larger than Sfe+i. 
After applying the projection map ws k+1> we have 

7T£ fc+1 (Reg afe=oc F) = ^^^(afc+i) , where 6 L(E k +i) ■ 

i 

The coefficients q t are functions of afc+2, . . . , a^, and can themselves be projected 
down to L(£fc + 2), • • ■ , £(£zv) in turn. Thus we can define: 
(44) RReg Qfc=oc F = tt Sjv o . . . o 7r Sfc+2 o 7r Sfc+1 (Reg Qfc=00 F) . 

Therefore RReg afc=00 F lies in C(a k+ 2, ■ ■ ■ ,ajv) <& L(R k+ i). 

Theorem 31. is unramified {definition ] 16]) , then the coefficients ofKKeg ak=00 F 
are multiple zeta values. 

The restricted regularization RReg a / for any element / G L(R k ) can be 
computed algorithmically by succesive differentation with respect to ak+i, ■ ■ ■ , c\n- 
This uses an induction on the weight which is illustrated on some examples in ^6.1.11 
below. The only difficulty in practice is to compute the constant terms 

R e S ajv =0 ° • • • ° R e ga fc + 1 =o(R e Sa fc = oo^) G C . 

The statement of the theorem is that this number lies in Z under suitable conditions. 
The proof uses associators in higher dimensions generalizing the ideas behind §5.1.31 

5.5. The integration algorithm. We wish to compute 

/>oo P oo 

1=1 ' ' ' / /i da\ . . .dctN , where f\ G R\ . 
Jo Jo 

The idea is to integrate, one variable at a time, using the poly logarithm functions 
in L(Ri). At the fc th stage of the integration process, we will have: 

/•OO poo 

1=1 ■ ■ ■ f k da k ■ ■ ■ da N , where f k G Z ® Q L(R k ) , 
Jo Jo 

and fk has weight at most k. The integrands fk are calculated recursively as follows: 

(1) The function f k has a primitive F k G Rk®R k L{Rk) of weight at most k + 1. 

(2) Since we can assume that Keg ak=0 Fk = 0, we have 

poo 

/ fkdak = Reg ak=00 F k . 
Jo 

In general, however, Heg ak=00 Fk is not an element of Z ®q L{R^i), so we 
must use restricted regularization instead. 

(3) Therefore, we define 

jfc+i = RReg afc=00 F fe G Z (g) Q L(Rk+i) . 
Although Jp 00 fk dak is not in general equal to fk+i, one can prove that 

/•oo poo poo poo 

/ • • ■ / fk da k ■ ■ ■ da N = I ■ I f k +i da k+ i . . ■ da N , 
Jo Jo Jo Jo 
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and the induction goes through. 
At the last stage of the integration, we have 

/>oo 

1=1 In da N = Reg aN=OQ F N , 
Jo 

where Fn is the primitive of Jn defined in (2) above. Since Fn S Z <E> L(Rn), its 
regularized value at oo is a Q-linear combination of multiple zeta values by £ )5.1.3I 
This proves that / S Z. 

Since the weight of the integrand increases by at most one at each integration, 
we can use the above argument to bound the weights of the periods we obtain. 

Remark 32. There is never any need to choose branches of any hyper logarithm 
functions, because one can prove that every function which occurs in the above 
process is in fact holomorphic on the open hypercube 

{(afc, . . . , a N ) 6 M. N ~ k : < a k , . . . ,a N < oo} . 

The functions fk can, however, have logarithmic singularities along the codimension 
1 faces of this hypercube. A detailed study of the locus of singularities in the 
Feynman case, and the corresponding compactifications, will be given in [9]. 

Note that the entire integration process is algorithmic, from the calculation of 
primitives to the regularizations at infinity, and can be reduced to a sequence of 
elementary manipulations on symbols. 
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6. Examples of hyperlogarithms and regularization 

We consider some one-dimensional examples which will occur in the calculation 
of the wheel with three spokes diagram. First let N = 1, gq = 0, and u\ — —1, and 
denote the corresponding alphabet by X = {xo,xi}. We set: 

(45) E a = {0, -1} , and = Q \x, -, — |— 

L x x + 1- 

and L(E X ) is the Os^. -module spanned by the hyperlogarithms ~Li w (x), where w is 
any word in X. Thus, in L(Y, X ), there are exactly two hyperlogarithms of weight 1: 

Li xo = log(ac) , Li xi = log (a; + 1) . 

Note the absence of a minus sign in front of Li Xl , which differs from the usual 
convention Lii(z) = — log(l — z). In weight two there are precisely four functions: 

Li x 2(x) ,Li X0Xl (a:) , Li XlXo (x) , Li x 2(x) . 

Using definition E01 one verifies that these can also be expressed as: 

ilog 2 (x) , -Li 2 (-x) , -Li 2 (x + 1) , ^\og 2 (x + l) , 

respectively. They are related by a single shuffle product: 

Li Xo (x)Li Xl (x) = Li XoXl (x) + Li XlXf) (x) . 

Finally in weight three there are exactly 8 linearly independent hyperlogarithms. 
One can check, for example, that Li x 2 Xi (x) = — Li3(— x). If we denote the regularized 
values at infinity by Coo(w) = ^ e Sx=oc^w(x), we have from j j5.1.3[ 

Coo(x^) = , Coo(x xi) - -C(2) , Coo(xiXo) - C(2) , Coo(x Xl) = . 

Likewise, in weight three, we deduce from lemma [27] that: 

(46) Cco(xix oXl ) = -2C(3) , Coo(x x?) = C(3) , Coo(x?x ) = C(3) 

Coo(XqXi) = ,Coo(x XlX ) = , Coo(xiX^) = , Coo(Xq) = , Coo(x?) = . 

6.0.1. The enlarged case {0, ±1}. Now consider the case when N — 2, o§ = 0, 
u\ = —1, 02 = 1. We set 

11 1 " 

x' X — 1 ' X + 1- 

We denote the corresponding alphabet by X + = {x ,xi,x 2 }. The algebra L(E+) 
now contains 3 hyperlogarithms in weight 1, namely: 

Li Xo (x) = log(x) , Li Xl (x) = log(x + 1) , Li X2 (x) = log(x - 1) , 

and 3™ hyperlogarithms in weight n. In this case, the regularized values at infinity 
are no longer multiple zeta values {Z 1 ), but the larger set of alternating multiple 
zeta values (Z 2 ). Via the inclusion L(E X ) — ► L(E+), we can identify L(E X ) as a 
sub-algebra of L(E+). In i )5.4l we defined a projection map: 

tt Sx : L(E+) — O s + £(E») 

f(x)L w (x) i > f(x)L w (x) ifw£l x , 

f(x)L w (x) i > if w contains the letter x 2 . 



(47) E+ = {0,-1,1}, and O s + = 



x, • 
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6.1. Second example. We now consider an example in two variables. Let 



(48) T, y = {<7o,<7i,cr 2 ,<7 3 } , where cr = , a x 
We set 

-- 



1 

y, -, 



i 



i 



■1 , Cr 2 
1 



-X , (T 3 



X+l 



y y+1 x + y y + ^ . 
Let us denote the corresponding alphabet by Y = {yo, yi, Y2, Y3}, and thus 
L yo(y) = l °g(y) , L yi (y) = log(y + 1) , L y . 2 (y) = log(x + y) - \og{x) , 

L vs (y) = l °s( x y + x + y)- \og{x) . 

Thus L(T, y ) is spanned (as a vector space) in weight n by exactly 4™ hyperlogarithm 
functions L w (y), which can be considered as functions of the single variable y, for 
constant x. If we wish to consider the dependence on x, we must remove the 
singular locus where the Oi collide. 

Recall that 0\ was defined to be 0£ B [(<7j), (<7j — <Jj) _1 ], and so 



(49) 



x x + l x — 1 



1 1 



y- 



y' y+ 1' x + y' y+ ^-J 



This is geometrically a linear fibration over O s + and is pictured below: 



y = oo 



y = 




Let [/ = C 2 \{a; = 0, x = ±1, y = 0, y — — 1, .x + y = 0, + a; + y = 0}, as 
shown. The elements of L(T, y ) can be viewed as multi-valued functions on U. The 
full space of polylogarithms is a product of two hyperlogarithm algebras: 

L(0+ ) = L(£+)L(£ y ) . 

It is generated by elements f(x,y) L Wl (x)L W2 (y) where w\ € X x , W2 <G Y x , and 
f(x, y) G 0^ . The set of elements of weight two are: 

Lixix^z) , Li Xi (x)Li ys (y) , Li yiyj (y) , 

where the x's are in X + , and the y's are in Y. There are, respectively, 9, 12, and 
16 such elements, giving a total of 37 polylogarithms of weight 2. 
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6.1.1. Regularized values. Let Z denote the Q-algebra spanned by all multiple zeta 
values. We defined a regularization map 

Reg y=00 : — > Z <8>q L(0£ e ) • 

The regularized values at y = of all functions L w (y), w GY X vanish. At y = oo, 

(50) Reg y=oa L yo {y) = , Reg y=oc L yi (y) = , 

^eg y=oc L y2 (y) = - log(x) , Reg y=oc L y3 (y) = log (a; + 1) - log(x) . 

We now compute the regularised values at infinity of some functions of weight 2. 
To compute Reg y=oc L w (y), we can differentiate with respect to x, which gives a 
function of lower weight, compute the regularised value at infinity of this function 
by induction, and take the primitive with respect to x. In other words, 

^g y=oc L w (y) = J Rcg y=00 {^-L w (y)^j dx . 

The constant of integration is determined from the regularized values: 

^eg y=oo Reg x=0 L w (y) 6 Z , 
which can be calculated by a generalization of the associator argument of M5.1.3I 

6.1.2. First example. Let w = yoY2- To compute -^L w (y), we have: 

d d d 8 5 11/1 In 

Q-y^yoyAv) - Q^Q- y L y^y) - q^-^Av) - -{—^ -) ■ 

Thus 



Tx L ^ M =l l^Ty) dy = --x LM 



The constant of integration is determined by the fact that L yoy2 (y) , and hence 
^L yo y 2 (y), vanishes along y = 0. Taking the regularized value at y — oo gives 

Reg y=00 (J^L yoy2 (y)) = -^Reg y=oc L y2 (y) = . 

By taking a primitive with respect to x (working now in we deduce that 

/log x 1 
—^-dx = L x 2 (x)+c=- \og{x) 2 + c , 

where c satisfies c = Reg x=Q Reg y=00 L yoy2 (y) = 0. It follows that 
(51) Reg a=oc i yoy2 (?/) = L x 2(x) . 

6.1.3. Second example. Let w = y3yi- By a similar calculation, we have 

-EZ L ysyi(y) = — V + 1 t L yi(v) TT L ys(v) ■ 

ox xy + x + y x + 1 

Taking the regularized limit at y = oo gives 

d 1 

^R-eg y=0O i y3yi (2/) = — y(log(a;) - log(ir + 1)) . 

To determine the constant of integration, we use the fact that Reg 2 , =0 Reg a=oo Z/ y3yi (y) 
C(2). Taking a primitive with respect to x in the ring L(E+), we deduce that 

Regv=oo i y 3 yi(j/) = L ^X X ) - L ^( x ) + C( 2 ) • 
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6.1.4. Further examples. In a similar way, one can verify the following: 
(52) Reg y=00 L yoy3 0) = L x 2 (x) + L x 2 (x) - L XoXl (x) - L XlXo (x) , 

Reg y=00 L y3y2 (y) = L XoXl (x) - L x 2 (x) + L x 2 (x) , 
which will be used in the calculation of the wheel with three spokes diagram overleaf. 

6.1.5. Restricted regularization. All the previous examples happen to lie in L(T, X ). 
In the general case, one obtains an answer with an extra singularity at x = 1, i.e., 
a hyperlogarithm in L(E+). For example, one can check that: 



( 53 ) in L M*{y) 



dx yiyAy) ~ x(x-l) x-1 



^g y=00 L yiy2 (y) = L X2Xo {x) . 

The reason this happens is because L yiy2 (y) has a singularity at o\ = a 2 - In this 
case, the restricted regularization map of §5.41 gives 

RReg y=oc i yiy2 (j/) = ttz x L X2Xo (x) = . 

Since the regularized value of every hyperlogarithm occurring in the calculation of 
the wheel with three spokes diagram overleaf already lies in L(T, X ), we do not in 
fact need to use the regularized restriction in the calculation. 

Remark 33. In order to compute the values Heg x=0 Reg y=OD L w (y), one can set 

t\ = -0-2 , h = -0-3 , and t 3 = -y . 

The hyperlogarithm L w (y), viewed as a function of the independent variables 
^1,^2,^3, is then a unipotent function on the moduli space 

m „ r 1111 1 1 1 1 1 

UH ,6 = bpecZ h,t 2 ,t 3 , 



L t\ t2 £3 1 — t\ 1 — t<z 1 — £3 t\ — t<z t\ — t 3 t2—t 3 - 

Setting t\ = x,t2 = j^j, ^3 = — V defines a surface inside 9Jto,6, an d taking the limit 
first as y — > and then x — » 0, corresponds to taking the regularized limit of L w (y) 
at the tangential base-point at t\ — 0, t% — 0, t 3 — defined by the sector 

< t 3 < t 2 < h < 1 . 

Equivalently, this corresponds to a unique point on the compactification 9Jto,6 which 
lies in the deepest stratum, and defined over the integers. Likewise, taking the 
regularized limit at x — 0, y — 00 corresponds to taking a limit at a different point 
in DJIq.q. By the results of [Sj, this limit can be described by a generalized associator 
and is expressible in terms of multiple zeta values. The general case is similar and 
follows from the ramification condition (|33[) (this will be discussed in [9]). 
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7. The wheel with three spokes 

7.1. The Fubini reduction algorithm. Let G be the wheel with three spokes 
depicted in fig. 1. We take A = 5, and set a$ = 1 from the outset. We therefore 
have S' — {Ug\ a5= i}, where U Q = Uccte + Va, and Ug, Vg are given by $8§. We 
wish to compute integrals of the form: 

j f ntiiog""(^)iog m (t/ 5 ) , . . , , 

(54) ±g = I —2 —aaedaida2da3aa4 , ra^m G N , 

Ja 5 =i <~>~ 

by integrating with respect to the variables ol\, a 2 , o 3 , 04, in that order. The 
Fubini reduction algorithm gives: 

Srgi = {oio 2 + O1O4 + a>ia 3 + o 3 o 2 + o 3 o 4 + o 2 + a 3 + o 4 , 

010304 + 010204 + 030204 + 0103 + 0203 + a 2 a^ + u\U/± + aia 2 } . 

After reducing with respect to Oi, we have: 

^[6,11 = { a2 + a 3 + 04 + 0203 + 0304 , 02 + «3 + 04 , 0304 + "3 + 04 , 

0204 + a3«4 + 2 + Qf3 + 04 , tt3tt4 + 2 + 3 + 04 } 

Reducing with respect to the variable a 2 gives: 

^[6,1,2] = i" 4 + 1 ' a 3 + 1 , "3«4 + o 3 + 0:4 , a 3 + o 4 } . 
Finally, by reducing with respect to the variable 03, we obtain: 

S [6A,2,3] = {"4 + 1} ■ 

Note that at each stage, every polynomial is linear with respect to every variable. 
The corresponding sets of singularities Si, ... ,£4 are therefore: 

trr\ r. r n a 2 + a 3 + 04 + a 2 a 3 + a 3 a 4 0203 + 0204 + 020304 

(55) Si = {0, , }, 

0:2+0:3+ 04 02 + «3 + 04 + 02«4 + O3O4 

^ r n / x / v 03+04+0304 03 + 04 + 0:304 

S 2 = {0, -(o 3 + o 4 ), -(o 3 + o 4 + a 3 a 4 ), — , — } , 

o 3 + 1 04 + 1 

S 3 = {0,-1, -04, , S 4 = {0,-1}. 

04 + 1 

Taking the limits as 02 — ► 0, 03 — > 0, 04 — > 0, in that order, we obtain singularities 
in 0, —1 only. Therefore the conditions of theorem [T8l are satisfied. 

Corollary 34. Every integral of the form \5Ji\ lies in Z. 

We have a nested sequence of rings Q C R4 C R3 C R2 C Ri, defined as follows: 

R 4 = Q[o4,a7 1 ,(a 4 + l)" 1 ] , 

R 3 = i? 4 [o 3 , 03 X , (a 3 + l) -1 , (a 3 + a 4 ) _1 , (030:4 + a 3 + o 4 ) _1 ] , 

R2 = i?3[o2, o^ 1 , (a 2 + o 3 + 0:4 + a 2 03 + a3a 4 ) _1 , (a 2 + a 3 + a 4 ) _1 , 

(02O4 + 0304 + o 2 + 03 + o 4 ) _1 , (03O4 + o 2 + O3 + 04) _1 ] , 
Rl = #2[oi, Oi , (0102 + 0104 + 0301 + O3Q2 + 030:4 + 0:2+03 + 4 ) _1 ] , 

(030104 + 030204 + 040:204 + a\a 3 + 0203 + 0204 + 0104 + 0102) ~ ] . 

One can compute any particular integral (|54[) by working in the algebra of polylog- 
arithms L(Ri) = L(Si) L(S 2 ) L(S 3 ) L(S 4 ), where L(Si), . . . , L(S 4 ) are hyperlog- 
arithm algebras on 3, 5, 4, 2 letters respectively. 
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7.2. Calculation of the leading term for the wheel with 3 spokes. We 

illustrate the method by calculating in complete detail the Feynman amplitude of 
the wheel with three spokes G, and reprove the result, due originally to Broadhurst 
and Kreimer, that it evaluates to 6 £(3). 

Let U = U G , V = Vg be as in ©. We set A = {5}, and wish to compute (123]) : 

daida 2 da 3 da 4 . 

{ a ,=i} uv 

From equation ([23)1 . we have D = a 2 a^ + a 4 «5 + a 3 a 4 + a 3 a 5 , and 

r/(l,2) _ a3Q , 4 _)_ a 3 a 5 + a 4 a 5 , = a 3 + a 5 , = a 3 + a 4 , U\2 = 1 . 

yf 1 ' 2 ) = , V% — — ct 3 a 4 + a 3 a 5 + a 4 a 5 , V\% = a 4 + a 5 . 
Let us set a$ — 1 once and for all. Using the notation ([26]) we have: 



rrrW rr(i,2)in n(2)\ (a 3 + 1) log(«3 + 1) . log(a 3 Q4 + t*3 + 0*4) 

1^2 > U \ u 2> u s - — 7 ; 1 H 1 ; ) 

a 3 {a 3 a 4 + a 3 + a 4 ) a 3 a 4 + a 3 + a 4 

fjj rr (2)i n n f2h log(a 3 a 4 + a 3 + a 4 ) - log(a 3 + a 4 ) log(a 3 + a 4 ) 
1^12, <A 1^2,^ } = ^7 1 

{V 2 W ,V^\D 2 ,D^} 



CX3O14 (a 3 a 4 + a 3 + a 4 ) 

2 log(a 3 a 4 + a 3 + a 4 ) 



(a 3 a 4 + a 3 + a 4 ) 

n/ t/(2)m n(2)i _ («4 + 1) log(a 4 + 1) log(a 3 a 4 + a 3 + a 4 ) 

{V12, Vi \U%,1J^ >) - — -, ; ; H ; ; • 

a 4 {a 3 a 4 + a 3 + a 4 ) a 3 c*4 + a 3 + a 4 

By corollary [51 we can skip the first two integration steps and go straight to 

(56) Ig = [ [ ( log(a 3 + 1) + log(a 3 + a 4 ) - log(a 3 a 4 + a 3 + a 4 )) 

Jo Jo a 3 a 4 V / 

H -. s(( a 4 + 1) log(a 4 + 1) — log(«3 + 1) — a 4 log(a 3 + a 4 ))da 3 da 4 • 

a 4 (a 3 a 4 + a 3 + a 4 ) V / 

One can verify that the integrand has no polar singularities along the faces of 
the domain of integration X = [0, 00] x [0, 00]. It has a pole of order exactly 
one along the hypersurface a 3 a 4 + a 3 + a 4 = 0, which meets the boundary of X in 
codimension two. Even though the integral converges, the simplest way to compute 
it is to calculate the regularised integral of each individual term, which may be at 
most logarithmically divergent, and add the contributions together. 

We set a 3 = y, a 4 = x. In R 3 and R 4 were called Os y and 0£ x , respectively. 
In the hyperlogarithm notation of tj5.ll the integrand of Iq is an element of L{R 4 ) = 
L{Oy, v ) of weight 1. We rewrite Iq as follows: 

!g= f f -j- (l Yi (y) + L Y2 (y) - L Y3 (y)) + -L^ (x) ( 1 - 
Jo Jo X V y / x 

7— TTV „ , x ( L yi (v) + xL Y2 (v) + xL *o 0)) dxdy . 

x(x + 1) y + ^ V / 

Using theorem 1301 we can take a primitive with respect to y, which gives: 

f°° 1 r 

Ig = / -Rcg y=oc L yoyi (y) + L yoY2 (y) - i yoy3 (y) + L Xl (x)L y3 (y) 
Jo x L 

~x{xTT) RcSv=00 Lysn ^ + xLy3Y2 ^ + xLx " ^ Ly3 ^ dx ' 
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Using the shuffle relations, and equation (1501) . observe that 
Reg y=oc L Xl {x)L y3 (y) = L Xl (x) (L Xl (x) - L Xo (x)) = 2 L x 2 (x) - L x 
and similarly, 

R-eg y=00 L xo (x)L y3 (y) = L XQX1 (ar) + L X1X0 (x) - 2 L x 2 (x) . 
All the remaining regularized limits were given in tj6.ll Substituting them in gives: 
" 1 



It 



G 



1 



C(2) + L XQXl (x) + L XlXo (x) - L x 2 (x) + L Xl (x) (L Xl (x) - L XQ (x)) 



C(2)+L XlXo (x)-L x 2(x)+x(L X(>Xl (x)-L x 2 (x)+L x 2 o (x)+L Xo (x)(L Xl (x)-L Xo (x))) 



x(x + 1) 

After decomposing into partial fractions, and expanding out the shuffle products, 
we obtain the fourth step: 



(57) I G = 



1 



2L x 2 i (x) - L XlX0 (x)) 



1 



1 



C(2) + L x? (x)-2L X0Xl (x) 



dx 



This is an integral of functions of weight at most two in a single variable x, in 
the hyperlogarithm algebra Z ®q L{R±) = Z ®q L(E x ). A priori, the regularized 
limits at y = oo which were substituted in at the previous stage could also have 
had singularities in £+, i.e., at x — 1 also, but the method (the Fubini argument) 
predicts that the integrand at the fourth stage (|57|) cannot, as is indeed the case. 
To complete the calculation, we work in L(Ri). A further integration gives: 

I G = Reg x=00 (2 L X()X 2 (x) - L XoXlXB (x) + ((2)L X1 (x) + L XlX 2 (x) - 2 i XlXoXl (xfj , 

Therefore at the fifth and final step, we obtain 

IG = 2Coc(x Xi) - Coo(x XlX ) + Coo(xiXo) - 2Coo(xiX Xi) . 

Substituting the values given in l|46p . we conclude that 

/g = 6£(3) . 



Remark 35. One can rewrite (|57[) using only dilogarithms and logarithms: 

(58) I G = [ — (-Li 2 (-a4)+log 2 (a4 + l)-log(a4)log(a4 + l)) 
Jo «4 \ / 

+ ^-(2Li 2 (-a 4 ) + Jlog 2 (a 4 ) +((2))da 4 , 

«4 + 1 V I J 

where each term in the integrand has singularities contained in {0, —1, 00}. One 
can obtain (|58p directly from (|56|) by integrating using the dilogarithm function, 
regularizing at infinity, and applying the inversion relation for Li 2 to arrive at (|58[) . 
Although this gives a substantial shortcut, such a method is ad hoc, and will not 
work in a more general setting. 
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